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Giant oxygen isotope effects observed in colossal magnetoresistance manganites are investigated 
by employing the combined model of the double exchange and interacting lattice polaron mechanism. 
We have shown that the isotope effects on Tc in the metallic phase and Tco in the charge ordered 
phase of manganites can be explained well in terms of the double exchange and polaron narrowing 
factors with reasonable physical parameters. 
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\^0 . A recent discovery of colossal magnetoresistance (CMR) phenomena has stimulated enormous interest in the per- 

ovskite manganites, R1_.zA3.MnO 3 (RAMO: R = rare-earth; A = divalent cation) [[ljfj]. They exhibit a very rich phase 
~ l , diagram depending on the doping concentration, temperature, and pressure: antiferromagnetic (AFM) insulator, fer- 
romagnetic (FM) metal, charge ordered (CO) insulator |3|. These novel features suggest that several interactions 
originating from the spin, charge, and lattice degrees of freedom are competing. For instance, the correlation between 
ferromagnetism and metallic conductivity for 0.2 < x < 0.5 was explained qualitatively in terms of the double exchange 
(DE) mechanism [Q. On the other hand, Jahn- Teller polaron effects due to a strong electron-phonon interaction are 
C$ ■ thought to be responsible for anomalous properties of manganites [pi . 

The most evident signature for the lattice degree of freedom is considered to be the isotope effect, because it can 
affect the electronic properties only through the electron-phonon interaction. Indeed many recent experiments, such as 
magnetization, resistivity, thermal expansion, electron paramagnetic resonance, and magnetostriction measurements, 
have revealed giant isotope effects both in the metallic and CO phases of RAMO p-fl-fl. The isotope effect on 
the magnetic phase transition temperature Tc in the metallic phase |3K| and on the charge ordering transition 
temperature Tco PHjjl] were confirmed by replacing 16 with 18 0. Zhao et al. [pi have obtained a shift of 20K in 
Tc in Lao.sCao^MnOa, yielding the isotope exponent cv c =0.85 (a c = — AlnTc/AlnMo, Mo- oxygen isotope mass). 
On the other hand, Franck et al. Q obtained a c ~ 0.4 for the same Ca concentration. This large difference in a c 
between two groups is suspected to arise from the different sample stoichiometry due to oxygen excess ||. 

For the CO phases, Isaac et al. |ll]] observed ATco with a negative oxygen isotope exponent a co = —0.41 
in Lao.57Cao.43Mn03. Further, novel crossovers from a metallic to a CO insulating ground state are observed in 
(Lao.5Ndo.5)2/3Cai/ 3 Mn03 |l(j,[l^] and in (Lao.i75Pr .525)Cao.3Mn03 Jl5| by the oxygen isotope exchange. The above 
systems are metallic with 16 0, but are very close to the phase boundary between the FM metal and the CO insulator. 
A more interesting feature is manifested in the CO phases of Ndo.sSro.sMnOa and Lao.sCao.sMnOs 13 13], which 



show the strong magnetic field dependence of the isotope effect. Under the magnetic field of 5.4 T, ATco increases 
substantially from 21K to 43K for Ndo.sSro.sMnOa, and from 9K to 40K for Lao.sCao.sMnOa, respectively. Thus, 

S the magnitude of a co (= — A In Tco/ A In Mo) increases rapidly with the magnetic field. This feature is in contrast 

to the case of a c in the metallic phase that is rather insensitive to the magnetic field J^,|ljl . It has been pointed out 
that these anomalous features are difficult to understand within the existing theories of the charge ordering transition 

c: mi. 

The isotope effects in the metallic phase can be explained properly by using the small polaron model M. In 
contrast, the understanding of the isotope effects in the CO phase is not straightforward. In the previous study p6[ , 
we have described successfully the lattice and magnetic properties of the CO phase of half-doped Mn-oxides using the 
combined model of the double exchange and interacting lattice polaron mechanism. The electron lattice of the CO 
state was viewed as the generalized Wigner crystal, that is, the CO state becomes stable when the repulsive Coulomb 
interaction between carriers dominates over the kinetic energy of carriers. In this paper, we have explored the isotope 
effects in the CO phase of CMR manganites on the basis of the similar model. We have found that the isotope effects 
and related features in the CO phase of manganites are well described in terms of the present model. 

Let us first examine the isotope effect in the metallic phase. The standard small polaron theory shows that the 
strong electron-phonon interaction reduces the electron hopping parameter by the polaron narrowing factor |17|] , 



t = t exp 



■5> q | 2 (A 9 + l/2) 



(1) 



with u q = M q /u q (exp(iq • 5) — 1). Here M q is the parameter of the electron-phonon interaction strength, uj q is the 



phonon frequency, N q is the boson distribution function, and S is the displacement to the neighboring site. The 
electron-phonon interaction parameter M q depends both on the ion mass and the phonon frequency in the form of 
M q ~ [MionOj)^ 1 / 2 . Hence ^ |w q | 2 is proportional to \JMi on with the assumption of uj ~ (Mi on ) -1 / 2 . In the DE 
model, Tc is proportional to the hopping parameter |18|, i.e., Tc ~ i ~ iexp [— jy/Mion] , and so the material with 



heavier ion mass would have the lower Tc- This is nothing but the isotope effect observed in the metallic phase of 
RAMO. Note that, with increasing ^ |w q | 2 , one gets the larger isotope effect. 

The experimental results of Franck et al. (a c — 0.34) J9| can be reproduced with ^ \uq} 2 — 1-2 and u> — 0.07eV. In 
the same framework, one can also explain the reduction of a c with increasing the tolerance factor P,p| . The increase of 
the tolerance factor leads to the increase of the Mn— O— Mn bond angle and the decrease of the bond length between 
neighboring Mn sites p9 |. One can thus expect that the phonon frequency becomes hardened due to reduced bond 
length (g0|, and accordingly ^ |u q | 2 decreases considering ^ |w q | 2 ~ (<^)~ 3 - Therefore, we have smaller a c with 
increasing the tolerance factor in agreement with the observation. In addition, the increase of the Mn— O— Mn bond 
angle and the decrease of bond length will certainly increase the hopping parameter and so weakens the polaron 



formation |21|. This effect plays a role of reducing ^ l^ql 2 and a c further. 

As for the giant isotope effects in manganites, Nagaev |J2^] has proposed mechanisms which are not related to the 
electron-phonon interaction. He claimed that a change in the oxygen isotope leads inevitably to a change in the 
sample stoichiometry. So thermodynamic equilibrium densities of oxygen vacancies or excess oxygen atoms depend 
on the mass of the oxygen atom and so corresponding carrier density has isotope dependence. In addition, the oxygen 
nonstoichiometry (vacancies) produces the isotope dependent pressure by causing a change in the volume. These 
mechanisms seem to explain observed isotope effects in Lao.sCao^MnOa which are very sample dependent (a c =0.85 
vs. a c ~ 0.4), as mentioned above. This scenario, however, works only for the low carrier density regime (x < 3/8). 
Moreover, more recent experiments have shown that there are intrinsic isotope effects which are not caused by any 
possible differences in the oxygen content |]|l(| . Also it is hard to explain with these mechanisms the isotope effects 
in the CO phase of manganites. 

Now let us consider the isotope effect in the CO phase. Since Tqo is not linearly proportional to the hopping 
parameter in this case, one has to find the relation between Tco and the hopping parameter. The extended polaron 
Hamiltonian for the half-doped manganites incorporating the electron-phonon interaction is expressed as 

ij ' <ij> q jq 

where tfj(= i(cos|), 9: angle between two spins on neighboring sites) is the DE hopping parameter and Vb is the 
Coulomb repulsion between two electrons of neighboring sites. Using the polaron canonical transformation, one can 
get 

H = J2 iiA c i + 2 V E ntTl 3 + E w « a q a q _ E n i A P' ( 3 ) 

ij <ij> q j 

where V ( = Vq — 2 ^2 —*• exp(iq • d) J corresponds to the renormalized Coulomb repulsion and A p (= ^2 — s -) is the 

polaron binding energy. The effective hopping t D is defined as in Eq. (fiT) with the bare hopping t replaced by t D . It 
can be shown that V does not depend on the ion mass when u> ~ (Afi orl ) -1 ' 2 . Thus the isotope effects will be realized 
only through i D , if they exist. 

To describe the CO state in the mean field approximation, the size of the unit cell is doubled by introducing formally 
two identical sublattices. This procedure automatically implies that the CO state can be only of a checkerboard type 
and certain restrictions are put on the lattice symmetry. Also, for simplicity, the electron hopping will be permitted 
only between sites belonging to different sublattices. Then the electronic part of the above Hamiltonian is reduced to 
the following form: 

H = E C4«/3 + 2W^(-ir4a a , (4) 

where a and (3 are sublattice indices, z is the number of nearest neighbors, and r is the charge ordering order parameter 
(—0.5 < t < 0.5), describing the deviation in electron concentration in the CO phase, i.e., n a = i + r; rip = i — r . 

Diagonalizing the above Hamiltonian, one gets two states with energy ±e(fc), where e(fc) = \/\i D (k)\ 2 + (zVt) 2 
p3[. Therefore, in the CO phase, the initial band is split into two bands with a charge gap of 2A = 2zVt. The order 
parameter r can be determined from the following BCS-like self-consistent equation, 



zV 1 ^-^ tanh(/3e(fc)/2) 
l -~ n^ 7{k) ' (5) 

where (3 = r-^f ■ The wave vector summation in Eq. rtq) can be replaced by the energy integration with the appropriate 
density of states (DOS). In the case of three-dimensional (3D) system, the DOS without the charge ordering /9o(e) 
can be chosen as a semielliptic form, po(e) — -^p\l ' B 2 — e 2 , where B(= zi ) is a half of the bandwidth without 

the charge ordering. Then the self-consistent equation at T — reads 1 = — b t T [K I «/ ,^\ 5 J — E I vpqn) ) > 

where b = -M?, and K(k) and E(k) are elliptic integral of the first and the second kind, respectively. 

By expanding the elliptic integrals into series in both the narrow and wide band width limits, the analytic form of 

the charge gap can be obtained: A « &. (l - \ (^) 2 ) for -^ < 1, and A = 45e~(^ +1 ) for ^ > 1. Figure [I] 
provides the overall behavior of Ao as a function of the band width. It is seen that Ao diminishes with increasing the 
hopping parameter, suggesting that the transition between the CO and the charge disordered state can be described 
by the competition between the hopping and the intersite Coulomb repulsion parameter. 

The self-consistent equation Eq. (pf) also determines the charge ordering temperature Tco- At Tco, the charge 
gap is closed in the spectrum and the order parameter r is equal to zero. In the narrow band limit, ^ < 1, the 

tanh function in Eq. (||) can be expanded into the series. Then Tco is given by ksTco ~ ^j- ( 1 — § iw) ) ■ Hence 
for the dispersionless band, it follows that ksTco ~ 4 — 2 > ' n accor dance with the previous work [jl6|, In the 
opposite case of the wide band limit, -y ^> 1, Tco can be obtained through the procedure analogous to the BCS 
theory [[24J . Careful retaining of all the terms that do not vanish under ^r-y — — ► 00 shows that the BCS result is 
obtained despite the nonconstant DOS of the present band, ksTco = — Ao ~ 0.567Ao, where 7 is the Euler constant 
(~ 0.5772). Thus, the ratio of -j~x — d° es n °t change significantly with varying the ratio of -M*. 

Overall dependence of Tco on -^ is presented in Fig. g. Note that Tco is zV/4 at B = and diminishes 
as B becomes larger. The observed decrease of Tco and the concomitant collapse of the CO state in half-doped 
manganites by applying the magnetic field or the pressure can be qualitatively understood in this context | J25| . That 
is, the magnetic field or the pressure increases the bandwidth to induce the instability of the CO state. 

Now, based on the above relation of Tco with the effective hopping parameter i D , one can account for the isotope 
effects observed in the CO phase. Since i D has ion mass dependence in the expression of the polaron narrowing factor, 
the replacement of 16 with heavier 18 gives rise to a smaller t D and so yields enhanced Tco- Let us estimate the 
reasonable physical parameters producing the experimental isotope exponent a co . For a given Tco — 150K, one 
can have t D = 0.05eV, J^ |« q | 2 = 1-2, and the screened Coulomb repulsion V = 0.018eV. The bare hopping 
parameter t D in the CO phase is thought to be smaller than that in the metallic phase because of the different spin 
ordering fluctuations. With these parameters, one could get the isotope exponent a co — —0.44 that is consistent with 

"Til. 
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experiments 

The isotope induced crossover from a metallic to a CO insulating ground state |I(J,|l2|,[l5| can be explained in the same 
way. The heavier ion mass of 18 reduces t, and so both Tc and the insulator to metal transition temperature Tmi 
decrease. In contrast, Tco increases with decreasing t, so that the CO insulating ground state becomes stabilized 
by the isotope exchange. One can also understand the anomalous isotope effects observed in the CO phases of 
Ndo.5Sr .5Mn03 and Lao.sCao.sMnOs under the strong magnetic field [|l3|,0. In the presence of the magnetic field, 
t D increases due to the DE term (cos |) and this will in turn reduce Tco- Under this circumstance, the isotope 
substitution of 18 reduces i D much more rapidly than the case without the magnetic field because of larger DE 
term, and so enhances the isotope exponent a co . This effect is clearly demonstrated in Fig. which plots — a co as 
a function of the bare DE hopping parameter t D . The magnitude of a co becomes larger with increasing t D , which 
corresponds to the case of applying the magnetic field. Direct comparison with the experiments may be possible, once 
the precise functional relation between t D and the magnetic field is known. The variation of —a co with t D and the 
corresponding physical parameters seem qualitatively consistent with observations 13 13, indicating that the present 
argument describes the main physics of isotope effects in the CO phase of RAMO. 

It is noteworthy that similar arguments were provided recently by Babushkina et al. p6j in analyzing their isotope 
induced metal-insulator transition in (Lao.i75Pro.525)Cao.3Mn03. In their model, the isotope effect is caused by the 
modification of the effective hopping integrals due to the change of interatomic distance by lattice vibrations. They 
obtained the relation between Too and the effective hopping parameter, and speculated that the large isotope effect 
arises from the hopping parameter that is very close to a certain critical value, t c = zV/2. In fact, this result 
corresponds to the present one in the narrow band limit 4cl. So their model is valid only in this limit. Further, 
they did not consider the polaron effect, so that the isotope mass dependence of the effective hopping integral is weak, 



-7=, as opposed to the present case ~ exp —\/M . Accordingly, their model is difficult to apply to the analysis of 



the anomalous isotope effects in the presence of the magnetic field |l3|;[l4| . 

To summarize, we have studied the isotope effects both in the metallic phase and the CO phase of RAMO, using 
the combined model of the double exchange and the interacting lattice polaron mechanism. We have investigated the 
charge ordering transition in the mean-field approximation, and obtained the relation between Tco and the effective 
hopping parameter. From this, we have shown that the various isotope effects including the magnetic field induced 
enhancement of a co in the CO phase of half-doped manganites can be explained in a natural way. 
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FIG. 1. The zero temperature charge gap Ao as a function of the half-bandwidth B (= zi D ) for the model 3D system. The 
values are scaled by zV where z is the number of nearest neighbors and V is the renormalized Coulomb interaction between 
carriers. 
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FIG. 2. The charge ordering transition temperature Too as a function of half-bandwidth B. 
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FIG. 3. The absolute value of the oxygen isotope exponent (— a co ) as a function of the bare DE hopping parameter t D . 
Calculations are performed for the parameters, ^ |w q | = 1.2, V — 0.018eV, Tco — 150K, ui = 0.07eV. 



